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Abstract 

We study the long time asymptotics of the relaxation dynamics of the totally 
asymmetric simple exclusion process on a ring. Evaluating the asymptotic amplitudes 
of the local currents by the algebraic Bethe ansatz method, we find the relaxation 
times starting from the step and alternating initial conditions are governed by dif- 
ferent eigenvalues of the Markov matrix. In both cases, the scaling exponents of the 
leading asymptotic amplitudes with respect to the total number of sites are found to 
be — 1. We also study the asymptotics of correlation functions such as the emptiness 
formation probability. 

PACS numbers: 02.30.1k, 02.50.Ey, 05.70.Ln 



1 Introduction 

The asymmetric simple exclusion process (ASEP) is one of the most fundamental models 
in nonequilibrium statistical mechanics [1-5]. The ASEP is a stochastic interacting parti- 
cle system consisting of biased random walkers obeying the exclusion principle, and have 
applications to biology [6], traffic flows [7,8] and quantum dots [9], for example. Like the 
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Ising model in equilibrium statistical mechanics, the ASEP is nowadays a paradigm in 
nonequilibrium statistical mechanics in the sense that many exact methods are amenable 
to extract various interesting nontrivial facts. For example, the matrix product representa- 
tion [10-15] of the steady state revealed boundary induced phase transitions [16]. For the 
last ten years, due to the development of the random matrix theory [17-24], the current 
fluctuations in the infinite system were shown to satisfy the Tracy- Widom distribution. 

The Bethe ansatz method, which was originated as a traditional method to analyze 
quantum integrable models such as the Heisenberg XXZ chain, can also be applied to 
analyze the ASEP. This comes from the fact that the Markov matrix describing the dy- 
namics of the ASEP is equivalent to the Hamiltonian of an integrable spin chain. Utilizing 
this fact, the relaxation times and spectral gaps were examined [25-33], and the exact 
expression for cumulants of currents and large deviation functions were obtained [34-38]. 

One of the latest developments of the Bethe ansatz approach is the study of the full 
relaxation dynamics of the totally asymmetric simple exclusion process (TASEP) on a ring 
by the algebraic Bethe ansatz method [39] . Previously, by noting the equivalence between 
the Markov matrix of the TASEP and the Hamiltonain of an integrable spin chain, the 
algebraic Bethe ansatz derivation of the Bethe ansatz equation and the construction of 
the determinant representation of the scalar product [40-42] was done. However, the 
power of the algebraic Bethe ansatz method was not fully utilized to study the dynamics 
of the TASEP, due to the difficulties of evaluating the multipoint form factors and the 
overlap between the initial state and the Bethe vector. Formulating the dynamics of the 
TASEP by evaluating the form factors and the overlap from the algebraic Bethe ansatz, 
we examined the full relaxtion dynamics of the local densities and currents, and found the 
scaling exponents of the asymptotic amplitudes starting from the step initial condition for 
example. The Monte Carlo method can be employed to study the full relaxation dynamics 
as well, but it is difficult to study the details of the dynamics, the long time asymptotics 
for example. 

In this paper, we study the long time behavior of the relaxation dynamics to the steady 
state by the algebraic Bethe ansatz method. We particularly focus on the local currents, 
and the two fundamental initial conditions: the step and the alternating initial conditions. 
The step initial condition is the case which the half of the system is consecutively occupied 
by the particles and the other half is empty at the initial time. The alternating initial 
condition is the case which all odd sites are occupied while all even sites are empty. We 
examine the asymptotic amplitudes of the local currents, and find that in contrast to 
the step initial condition, the asymptotic amplitudes corresponding to the lowest excited 
states of the Markov matrix vanish for the alternating initial condition. In other words, the 
relaxation time of the alternating initial condition is not governed by the lowest excited 
states. Instead, the second lowest excited state determines the relaxation time. Our 
discovery suggests that the naive guess of the relaxation time by considering only the 
Markov matrix may lead to an incorrect result. 

This paper is organized as follows. In the next section, we review the basics of the 
algebraic Bethe ansatz and the scalar products. In section 3, we formulate the dynamics 
of the local densities, currents and the emptiness formation probability of the TASEP by 
the algebraic Bethe ansatz method. This is achieved by evaluating the form factors and 
the overlap between the initial state and arbitrary Bethe vector for the case of step and 
alternating initial conditions. The low excited states of the Markov matrix are described in 
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section 4. The long time asymptotics for the step initial condition is analyzed in section 5. 
In section 6, we analyze the alternating initial condition and extract interesting difference 
from the step initial condition. Section 5 is devoted to the conclusion. 

2 Algebraic Bethe ansatz of the TASEP 

In this and the next sections, we formulate the dynamics of the TASEP on a periodic ring 
by the algebraic Bethe ansatz method. In this section, we review the relation between the 
TASEP and the integrable spin chain, and the basics of the algebraic Bethe ansatz and 
the scalar products. 

2.1 The definition of the TASEP 

We consider the TASEP on a periodic ring with M sites and N particles. Since the 
particles obey the exclusion rule, each site can be occupied by at most one particle. The 
dynamical rule of the TASEP is as follows: during the time interval dt, a particle at a 
site j jumps to (j + l)th site with probability dt, if the (j + l)th site is vacant. For 
convenience, we associate a Boolean variable T{ to every site i to indicate whether a 
particle is present (r^ = 1) or not (t^ = 0). The probability of being in the (normalized) 
state \ti, . . . ,tm) is denoted as Pt(r\, . . . ,tm)- The time evolution of the state vector 
\ip(t)} = J2 Ti =o l Ptiji-, ■ ■ ■ , tm)\t\i ■ ■ ■ , tm) is subject to the master equation 

±\m)=M\m)- a) 

Here the Markov matrix A4 of the TASEP is defined by 

• M = E{<^i + l(^I + i- 1 )}' ( 2 ) 

where := (<rj ± ia y -)/2 and crj' y ' z are the Pauli matrices acting on the jth site. Here 
we interpret the occupied (r, = 1) and unoccupied (r, = 0) state with spin down (| ^ 
)) and up state (| ti))i respectively. For example, we interpret that <r||n, . . . ,tm) = 
(— l) Tj |n, . . . , tm). We denote the vacuum state (state with no particle) := |0, . . . , 0). 
Starting from any initial condition, the state of the TASEP converges to the steady state 
|<Sjv) ( n °t normalized) 

\Sn) = ^2 a mi a m 2 - ■ ■ a m N \ty, (3) 

l<mi<-"<mjv<M 

in the long time limit. The steady state is an eigenvector of the Markov matrix with zero 
eigenvalue 

M\S N ) = 0. (4) 
We also define the dual vacuum state (0| := (0, . . . , 0| and the left steady state vector 

{Sn\= ( n \ a mi a m2 •••< J m JV > ( 5 ) 

l<mi<---<mjv<M 
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which is also an eigenvector of the Markov matrix with zero eigenvalue 

{S N \M = 0. (6) 
The inner product between \Sn) and {Sn\ can be easily calculated as 

Z N := (S N \S N ) = N]{M _ Ny - (7) 
2.2 Algebraic Bethe ansatz 

From the spin chain point of view, the Markov matrix (2), which describes the dynamics of 
the TASEP, is exactly a Hamiltonian of a one-dimensional integrable quantum spin chain. 
Therefore one can use the exact methods to examine the dynamics of the TASEP. The 
Bethe ansatz is one of the most traditional methods to treat quantum integrable models. 
The algebraic Bethe ansatz [40-42] is one of the variants of the Bethe ansatz, which can 
construct the eigenvectors as well as the eigenvalues of the Markov matrix (Hamiltonian). 
Moreover, the algebraic Bethe ansatz enables us to calculate form factors, which are the 
basic ingredients of the physical quantities. Therefore we formulate the TASEP by the 
algebraic Bethe ansatz method. 

What plays a fundamental role in an integrable model is the L-operator acting on the 
jth site 

L(j\u) = ussj + n(ul — u~ Sj) + a~ ~cr+ + c + cr~, (8) 

where Sj = (l + cr|)/2 and rij = (1 — <tJ)/2 are the projection operator onto the empty and 
filled states at jth site, respectively (see figure 1 for a pictorial description). The symbols 
s and n without subscript mean they act on the auxiliary space. 
The L-operator satisfies the RLL relation 

R(u, v)(L(n\u) <g> L{n\v)) = (L(n\v) <g> L(n\u))R(u, v), (9) 

where R(u, v) is the i?-matrix 



R(u, v) 



(10) 



//(«,«) \ 

o g(v,u) l o 

g(v,u) 

V f(v,u) J 

f(v,u)=-^ — o' 9(v,u) = |~ 2 ' ( U ) 

U — V u — V 

From the RLL relation (9), it follows that the monodromy matrix which is defined as a 
product of L operators 

" A(u) B{u) 
C{u) D{u) 



satisfies the RTT relation 

R(u, v)(T(u) (8) T(v)) = (T(v) ® T(u))R{u, v). (13) 
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Figure 1: The elements of the L-operator (8). The L-operator is represented as two 
crossing arrows. The left and the up arrow represents an auxiliary and a quantum space 
respectively. The spins on the left and the right around a vertex denote the input and the 
output of the auxiliary space, and the ones on the bottom and the top denote the input and 
the output of the quantum space, respectively. For example, (t |(| |jL(j|n)| \) \ \)j = 1. 



From the intertwining relation (13), one immediately finds the transfer matrix 

T (u) = U - M trT(u) = u- M {A{u) + D(u)), (14) 
forms a commuting family 

[t(u),t(v)}=0. (15) 
The Markov matrix (2) is constructed from the transfer matrix (14) as 

(16) 

The elements of the RTT relation (13) give the commutation relations between the ele- 



M = k'wsrwi— 
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ments of the transfer matrix A(u), B(u), C(u) and D{u). Some of them are listed as 

C(u)B(v) = g(u, v){A{u)D{v) - A(v)D(u)}, (17) 

A(u)B(v) = f(u, v)B(v)A{u) + g(v, u)B(u)A(v), (18) 

D(u)B(v) = f(y, u)B(v)D{u) + g(u, v)B(u)D(v), (19) 

[B(u),B(v)] = [C(u),C(v)]=0. (20) 

The arbitrary TV-particle state \il>({ u })) (resp. its dual (iIj({u})\) (not normalized) with ./V 
spectral parameters {u} = {m, U2, ■ ■ ■ , un} is constructed by a multiple action of B (resp. 
C) operator on the vacuum state |0) (respectively, (f2|) 

N N 

mu})) = HB( Uj )\n), (^({«})| = <nin^)- ( 21 ) 

i=i j=i 
Utilizing (18), (19), (20) and the action of A(u) and D(u) on the vacuum state 

A( u )\n) = u M \n), D{u)\n) = ( u - u- 1 )™^}, (22) 

one can show that the state \ip({u})) is an eigenstate of the transfer matrix (14) 

t(vM({u})} =Q N (v,{u})mu})}, (23) 

N U 2 N 2 

&n(v, m) = n + (i - v ~ 2 ) M n ~2~ — 2 ' ( 24 ) 

3 =X 3 j=l 3 

if the spectral parameters {u} = {u\,U2, . . . , itjv} satisfy the Bethe ansatz equation 

JV 

^-u- k 2 )- M uf N = (-l) N ^X[uj\ (25) 

for fe = 1, 2, • • • , JV. One can also show 

mu})\r(v) = mu})\Q N (v,{n}), (26) 

under the constraint (25). The eigenvalue of the Markov matrix (2) is given by the spectral 
parameters as 

1 8 N 1 

M(W) = ^(1, {u})-Q N (v, {u})\ v=1 = ^T37- ( 27 ) 

j=i j 

The steady state |SV) ((«SW|) corresponds to the eigenstate with zero eigenvalue which is 
given by setting the spectral parameters as u\ = U2 = ■ ■ ■ = un = oo. 

The scalar product, which plays an important role in this paper, has the following 
determinant form [42] 

(«i«=( n n ^>} det ^ (28) 

l N>j>k>l k 3 N>l>n>l 1 n ) 
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where P is an N x N matrix with matrix elements 

{/ \ JV-l / \ -N+l- 



Vi \ v. 



(29) 



Taking {v} = {u}, one gets the determinant representation of the norm 

N N 



(^(W)I^(W)) = f[uf N+M - 1) n ^~^det N Q, (30) 

j=l l,n=l 1 ™ 

with N x N matrix 

JV-1 + (M-JV + l)u7 2 

Q 3 k = ~ 2 —S jk - (1 - S jh ). (31) 

1- Uj 

Note that by use of Sylvester's determinant theorem, one can reduce the determinant in 
the above to 



N N + [M _ N)uJ 2 N 1 _-2 

det^Q = rr — - — 2 Lj -\ l -Y — ; — s 2 ■ 32 

f = \ l-uf y j^N + {M-N)uf] 
3 Form factors and Overlap 

We first formulate the relaxation dynamics of the TASEP to see what we should evaluate. 
Then we evaluate the form factors and the overlap between the initial state and arbitrary 
Bethe vector in the determinant and factorized forms. We consider two fundamental initial 
conditions: the step and alternating initial conditions. 

3.1 Formulation of the relaxation dynamics 

We formulate the relaxation dynamics of the TASEP by the algebraic Bethe ansatz. The 
time evolution of the expectation value for the physical quantity A starting from an initial 
state | In) is defined as 

(A) t = (S N \Ae Mt \I N ), (33) 

where (SV| is the left steady state vector (5). This definition comes from the fact that 
the TASEP is a stochastic process, and the coefficient Pt(ri, . . . ,tm) of the state vector 
\ip{t)) = er"\Ifj) directly gives the probability of being in the state \t±, . . . , tm)- We 
decompose the quantity (33) by inserting the resolution of identity 

T = \S N ){S N \ sr^ I^H^al 
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Here a labels arbitrary eigenstates except for the steady state. Then we find the local 
densities (ni)t = (1 — Sj)t and currents (ji)t = ((1 — Si)si+i)t are respectively given by 

, v N \- e M ^((S N \^ a ) - (S N \ Si \^ a ))(i; a \I N } 
\ n i)t — 77 + / , 7; 17-7 i 

_ N{M - N) ^ e^° t «5 jy |s i+ i|V; a ) - (Sjvl^+ilVQXVg/jv) 

" m(m - 1) + ^ <V*> ' ( } 

In the same way, one can also make the same decomposition for the emptiness formation 
probability (EFP) EFP(i, k)t = (sjSj+i • • • Sj+fe_i)t, which gives the probability that the 
sequence of k sites (i, (i + 1), . . . , (i + k — l)th sites) are all unoccupied 

EFP(i,k) t 

_ (M - N)\(M - k)\ sr e M « t (S N \s i s i+l ---s l+k - l \^ a ){ilj a \I N ) 

M\(M-k-N)\ + ^ {rPM ■ 1 j 

From the expressions of the physical quantities (35), (36) and (37), one notices that what 
we should evaluate is the norm (ip a \ip a ), the form factors (SVIsiSi+i • • • Sj+fc-ilVO an d the 
overlap between the initial state and arbitrary Bethe vector (ip a \I^f). The norm can be 
obtained as a limit of the scalar product (28) in the determinant form (30). What remains 
to be done is the evaluation of the form factors and the overlap. 

3.2 Form factors 

The form factor for the local operators Sj ■ ■ ■ Si+k-i is explicitly given by the following 
determinant form 

(S N \si ■ ■ ■ s i+k -i\ip({u})) 

N N 

=n( i -»j 2 ) i+, - i n»f +i n srb**^"-". <»> 

j=l j=l N>l>n>l 1 n 

where the N x N matrix V is written as 

(M-k) _ S^ ( _ 1)n (M-k)l 2(j-l-n) (3g) 

V ol "2^ L ) n\(M-k-n)\ Ul ' [M) 

n=0 

for 1 < j < N - 1 and 

y(M-k) _ _ V (M-k)\ 2(n-7V+l) m 

Vm ~ L> 1 j n\(M-k-ny. Ul ' 1 Uj 

n=N-l y ' 

Note that the overlap between the steady state and the Bethe vector (Sn\^) is obtained 
by setting i = 1, k = in (38). We show (38) by applying the approach of [42]. First, let 
us denote the monodromy matrix constructed from the (M — k + l)-th, • • • , M-th sites 
and its matrix elements as 

w - n *c*o = ( c k M ) • (41) 
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By definition, one has 



A.M-k+\{u) B M -k+i(u) 
C M -k+i{u) D M -k+i(u) 

A M _ k (u) B M _ k (u) \ / us k 

C M _ k {u) D M _ k (u) ) I at ul-u" x s k 



(42) 



from which we get 



B M -k+i(u) = A M - k (u)a k + B M -k(u)(uI ~ u 1 s k ), (43) 
C M -k+i{u) = C M -k{u)us k + D M _ k (u)a+. (44) 

Acting both sides of (43) by s k from the left, and (44) from the right, one has 

s k B M ~k+i(u) = {u- u~ 1 )B AI _ k (u)s k , (45) 
C M -k+i(u)s k = us k C M -k{u). (46) 

Utilizing (45) and (46), one can calculate the following generalized form factor as 

m{v})\s lS2 --- Sk \iJ({u})) 

N-l N 

= (n\ J] C(« J -)C , (« JV )a?a 2 ---a fc fl(ui)JJfl(u i )|n> 

3=1 3=2 
N-l N 

= (n\ H C(v ] )C(v N ) SlS2 ---s kSl B(u 1 ) ~[[B( Uj )\n) 

3=1 3=2 

N-l N 

=(Q\ ] [ C(vj)v N siC M -i(v N )s 2 ■■■s k (ui- u^ 1 )B M -i(ui)s 1 B(uj)\Q) 

3=1 3=2 
N N N 

= J\vj(uj - uj 1 )^ Y[ CM-i(vj)s 2 s k JJ B M -i(uj)\tt) 

3=1 3=1 3=1 



N N N 

= n - ^ Ti ) fc ^i n c M-k(vj) n B M _ k ( Uj )\n). (47) 

3=1 3=1 3=1 

Utilizing the cyclic property s k = r fc_J SjT^~ k ,r = t(1) and the action of r on the Bethe 
vector, one gets 

(ip({v})\siS i+1 ■ ■ ■ s i+fe _i|V>(M)) 

= ^(W)|(T J - 1 Sl r 1 - 4 )(r J - 1 S2 r 1 - J ) ■ ■ ■ (t^t^I^M)) 

= ^(W)|r t - 1 SlS2 ... Sfc r 1 - J ^(W)) 

-2\ *-l JV JV JV 

n v H u J - ^ ri ) fc (^i n cw^o n £w-*(«i)i«>- (48) 



iz5 

l-v7 2 



3=1 3=1 3=1 
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The form factor (38) can be obtained by taking a limit of the generalized form factor (48). 
First, note that the steady state can be obtained as 



N 



\S N )= lim ~[B M ( Uj )\n), 

{u}— i-oo 

3=1 

AT 

(s N \= lim (n|TT(7 M («,■), (49) 

where Bm(u) = u~( M_1 ) Bm{u) and Cm{u) = u~^ m ~^Cm(u)- We rewrite the generalized 
form factor (48) as 

AT N 

{Q\ Y[ C M {vj)siS i+ i ■ ■ ■ s i+k _i 

3=1 3=1 

{ I — u~ 2 \ i ~ 1 N N N 

= Hi 1 " u 7 2 ) k ( n \ II CM- k (v 3 ) n BM-k(uj)\n). (so) 

V V i / 3=1 3=1 3=1 

Taking {v} — > oo, one can show [42] 

N N 

(S N \Y[B M ( Uj )\n) = l[ U l H -^—^det N V^ M \ (51) 



3=1 k=l N>l>n>l U l Un 



where 



n=0 



nUM - n)\ 

n=N-l v ' 

Taking the limit {v} — >■ oo in (50) and inserting (51) to the right hand side, we have 

N 



(S N \siS i+1 ■ ■ ■ Si +k ^iY\_B M (ui)\Q) 



N N 



j=l j=l N>l>n>l I n 

Changing from Bm(u) to Bm(u), one obtains the determinant representation for form 
factors (38). 

3.3 Overlap: Step initial condition 

What remains to be done is to evaluate the overlap (tp a \lN) between the initial state and 
arbitrary Bethe vector. First we consider the step initial condition (see figure 2 (a)) 
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(a) 



(b) 




Figure 2: (a) The step initial condition and (b) the alternating initial condition of the 
TASEP on a ring. 

where the half of the system is consecutively occupied by the particles and the other half 
is empty. By graphical description of the L-operator, it is easy to see (figure 3) that the 
(normalized) initial state |/jv) = 1 1, - - - , 1, 0, . . . , 0) is given by |/jv) = ^(l)^!^) ( see [43] 

N M-N 

for the XXZ spin chain). Note that this initial state is not the eigenstate of the Markov 
matrix. 

Then we find the overlap (ip({v})\I N ) = (^{{v})\B{l) N = (Y>(M)|^({1})) can be 
obtained as a limit {u} — > 1 of the scalar product formula (28). One finds 

/_i\JV N (v - v _1 ) M 1 

ww)i^)= i iv(ivWTT i- 2 v n 7i—2 

2 2 3=1 V j N>j>k>l U 3 U k 

x detN ( 1 1 1 \k + T\ 1 \k ] ■ ( 55 ) 
We can furthermore simplify the determinant as 



fl 1 \ ( 2 k\ 

detN 0^ + (TT^ =detN 0^ ^ Wf—lY v i 



=det 



2 k 



k 



N 



(1 - v]) k 



2 JV(JV+l)/2 



det N ((l-vj) N ~ k ) 



lEU(i-v 

^ N+1)l2 Y\ N>]> k>M-vl) 

njli(i-«?)" 



(56) 



We made column operation in the determinant in the second equality, and used the formula 
for the Vandermonde determinant 

det N (xf- k )= H (x k - Xj ), (57) 
N>j>k>\ 
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Figure 3: Graphical description of \In) = B(l) for the step initial condition. The 
spins on the bottom line represents the vacuum, and each row corresponds to the B 
operator. Setting the spectral parameter on the bottom line to be one, one finds the two 
southwest internal spins freeze since (t \(X \\L(l\u = 1)| \)\ t)i = 1 and (4- |(t 1 1-^(1 |w = 
1)1 4)1 t)i = 0- Repeating the same game makes one see that all spins freeze. The spins 
on the top line is the resultant state of the action of B(1) N on the vacuum. 

in the last equality. Inserting (56) into (55), one gets the following simple form for the 
case of the step initial condition 

N 

m{v})\I N ) = life - vj l ) M - N vf~ l . (58) 
i=i 

Note that this form holds for arbitrary filling. 



3.4 Overlap: Alternating initial condition 

We now evaluate the overlap for the case of the alternating initial condition \In) = 
1 1, 0, 1, 0, . . . ,1,0) (figure 2 (b)), which all odd sites are occupied while all even sites are 
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empty. We consider the half-filling case. We find the following simple form 

N N 

wm)\in) = - vj 1 ) n ( v H - !)• ( 59 ) 

j=l j,k=l 

Let us show (59). For convenience, we use the spin notation (spin up and down states 
correspond to occupied and empty sites respectively). First, by representing the left hand 
side of (59) graphically (figure 4), one finds 

N N 

{ijj({v})\ hhUU ■■■ Im-i^m) = Y[v j Y[(v j -vJ 1 )V n {v 1 ,v 2 ,--- ,vn), (60) 

3=2 j=l 

V N (v!,V2, ■■■ ,v N ) =(n\C{v N ) ■ ■ ■ C(v 2 )D( Vl )\ hUUU ■ ■ ■ ±m-i). (61) 

We focus on T>n(vi,v 2 , ■ ■ ■ ,vn). Again, by graphical representation (figure 5), we note 
the following recursive relation 

N N 

V N (±l,v 2 , ■ ■ ■ ,v N ) = Y[ v 3 T[( v 3 ~ v j 1 )' D N-i{v2- l ■■■ ,v N ). (62) 

3=3 3=2 

Since (^({^1)1 -UT24-3T4 ' ' ' \-M-i\m) is symmetric with respect to v\,v 2 , • ■ • ,vn 
([B(vi), B(vj)\ = 0), T>n{v\,v 2 , • • • , vn) must be of the form 

N 

V N (vi,v 2 , ■■■ ,v N ) = ~[\_ v J 1j7 n(v 1 ,v 2 , ■ ■ ■ ,v N ), (63) 

3=2 

where J f n(vi,v 2 , • • • , vn) is a symmetric polynomial of v±,v 2 , ■ ■ ■ , vjy. Utilizing (63), the 
recursive relation for T>n(vi,v 2 , ■ ■ ■ , vn) (62) becomes the one for Tn(vi, v 2 , ■ ■ ■ , vjy) 

N 

T N (±l,v 2 , ■ ■ ■ ,v N ) = JJ(u| - l)Jjv-i(w2, • • • ,vn)- (64) 

3=2 

By symmetry, (64) extends to 

N 

Fn{vi,v 2 , ■ ■ ■ ,v N )\ Vk=±1 = J\(v] - l)T N -i(vi, - ■ ■ ,v k ,--- ,v N ). (65) 

One finds 

N 

J 7 n(vi,v 2 , ■■■ ,v N ) = [ (vjvl - 1), (66) 

j,k=l 
j<k 

solves the recursive relation (65). Combining (60), (63) and (66), one gets the factor- 
ized polynomial expression for the overlap (59) between the alternating initial state and 
arbitrary Bethe vector. 
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Figure 4: Graphical description of (60) . 

4 Excited states 

In section 2 and 3, we have formulated the dynamics of the TASEP by the algebraic Bethe 
ansatz by evaluating ingredients such as the form factors and overlap. In this section, we 
review the results of the excitation spectrum of the TASEP on a ring [28, 29, 31]. 



4.1 Algorithm 

We review the algorithm of computing the excitation spectrum of the TASEP on a ring. 
To this end, we make change of variables of the spectral parameters from Uj to Zj = 
(Zj + l)/(Zj — 1)) in this section. The Bethe ansatz equation can be rewritten as 

N 



(1 



z2\N 



n 2 



i 



z 4 + i 



for k 



1,2, . . . , N . The eigenvalue of the Markov matrix becomes 

N 7 1 
M({Z}) = ^^J 1 - 

3=1 



(67) 



(68) 
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Figure 5: Graphical description of (62) . 

A simple algorithm was proposed [28] to calculate the roots of the Bethe ansatz equation. 
Noting the right hand side of (67) does not depend on the index k, we define a parameter 
u as 

(1 - Zl) N = -e™. (69) 

The roots of this equation are 



Z m = —ZN+m = \A — Urn, (70) 
y m = e 2w(u+i)/M e 4ni(m-l)/M ^ /y-g 

for m = 1, 2, • • ■ ,N, and {y} = {y 1 , y 2 , . . . , y^v} satisfy 

< arg(yi) < arg(y 2 ) < ■ ■ ■ arg(?Mr) < 2vr. (72) 

It was proposed in [28] that choosing a sequence of quantum numbers {c(l), c(2), . . . , c(N)} 
satisfying 1 < c(l) < c(2) • • • < c(N) < M and determining the parameter it from 

^ Z - 1 

e™ = 2 M TT y c(j) (73) 

jJi Mj) + 1 

self-consistently by numerical iteration, one gets the Bethe roots {Z c n\, Z c f 2 \, . . . , Z c ^}. 
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4.2 Excited states 



By numerical calculation, we observe the following types of Bethe roots are the three 
lowest excited states. We mean a lower- lying excited state by a state whose real part 
of its corresponding eigenvalue of the Markov matrix is closer to zero. Here, we impose 
the ansatz that the eigenvalues of the low-lying excited states behave as ln(ReTW) = 
a — /3lnM,P = 3/2 by the following reasons. Since the lowest excited states believed to 
be true behave in this way, the exponents (3 for the other excited states should not be 
/3 > 3/2 for M large enough since no crossing across the lowest excited states is allowed. 
Next, by estimating several low- lying excited states by numerical observations for M ~ 20 
and making finite size scaling analysis of them for large M, we find the exponents for all 
of those states are close to /3 = 3/2, which is the reason why we impose the ansatz. 

(1) Type I [28,29,31] 

The Bethe roots corresponding to the quantum numbers 

{c(j) = j for j = 1, ■ ■ ■ , JV - 1, c(N) = N + 1}, (74) 
{c(j) = j + 1 for j = 1, • • • , N - 1, c(N) = 2N}, (75) 

give the lowest excited states. The simplest fitting from M = 256, 512, 1024 gives ln(— M\ s t) = 
1.89793 - 1.503511nM, which implies the KPZ scaling, i.e. M = CM~ 3 / 2 for M » 1. 

(2) Type II 

The second lowest excited state is given by the quantum number 

{c(j) = j + 1 for j = 1, • • • , N - 2, c(N - 1) = N + 1, c(N) = 2N}, (76) 

for M large enough. Finite size scaling from M = 256, 512, 1024 shows ln(— A^nd) = 
2.81592 — 1.5057681nAf , which shows the KPZ scaling again. This state will be important 
for the case of the alternating initial condition. 

(3) Type III 

The Bethe roots associated with the following four quantum numbers 

{c(j) = j for j = 1, ■ ■ ■ , N - 1, c(N) = N + 2}, (77) 

{c(j) = j for j = 1, • • • ,N - 2, c(N -1) = N, c(N) = N + 1}, (78) 

{c(l) = 1, c(j) = j + 1 for j = 2, • • • , N - 1, c(N) = 2N}, (79) 

{c(j) =j + l for j = I,-- - , N- 1, c(N) = 2N- 1}, (80) 

give the third lowest excited states. Conducting the finite size scaling from M = 256, 512, 1024 
shows ln(-Re7W 3rd ) = 2.86964 - 1.503471nM, which also belongs to the KPZ scaling. 



5 Step initial condition 

In this and the next section, we examine the relaxation times by studying long time 
asymptotics of the local current and emptiness formation probability. We consider the 
step initial condition in this section. 

We first analyze the local current. In the long time, the local current behaves as 

{jk)t ^ M(M +A ^> MlStt 

+ A 2nd (j k )e M ^ t + A 3rd (j fe )e Re ^ i cos(Im.M3 rd t + S k ) as t -> oo. (81) 
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Table 1: Table of the asymptotic amplitudes vs total number of sites (step initial condi- 



tion). 


(i,k) 


\n[—A(EFP(i, k))] 




1.93798 - 1.493321nM 


(iV -1,2) 


1.35223 - 1.055921nM 


(N — 2, 3) 


1.42851 - 1.016831nM 


(JV-3,4) 


1.22155 - 0.991481nM 


(iV-4,5) 


0.84935 - 0.967721nM 


(N — 5,6) 


0.36629 - 0.942901nM 


(N - 6, 7) 


-0.19885 - 0.916131nM 


(N-7,8) 


-0.82903 - 0.887051nM 


(iV-8,9) 


-1.51332 -0.855491nM 


(AT — 9, 10) 


-2.24430 - 0.821361nM 



The simplest fitting from M = 256, 512, 1024 shows ln[A lst (j N )} = 0.854409-0. 99685791nM, 
ln[-A 2ll d(jiv)] = 1-1875099 - 0.98389211nM and ln[|A 3rd (jjv)|] = 2.20350 - 0.983451nM 
from which one concludes ^4i s t(j7v) oc M , A 2n d(j7v) oc M~ l and | ^3rd ( J at ) I oc M~ l . 
Especially, Ai st (jjv) oc M _1 means that A\ st {jN) ^ as long as the total number of 
sites M is finite, confirming the relaxation time of the step initial condition r s t ep is de- 
termined by the lowest eigenvalue of the Markov matrix which shows the KPZ scaling 
Tstep — -Re(l/7Wi st ) oc M 3 / 2 . 

We can also examine the asymptotic amplitudes of the EFP. 

EFP{h k)t -+ MKM-fc M AT^f + A ( EFP h k » eMlStt as * -+ °°- ( 82 ) 

Table 2 is the results of the finite size scaling of the asymptotic amplitudes. One observes 
A(EFP(i,k)) oc M~ ak , and becomes smaller as the length k becomes longer. 

6 Alternating initial condition 

Next we study the alternating initial condition. Surprisingly, we find that the asymptotic 
amplitudes of the local currents associated with the lowest (Type I) and the third lowest 
(Type III) excited states vanish. This can be shown as follows. Rewriting the overlap (59) 
in terms of the spectral parameters Z c ^ in section 4, the following terms 

JV 

[] (Z c{j) +Z c{k) ), (83) 

j<k 

appear. Let us take a look at one of the excited states of Type I {c(j) = j for j = 1, • • • , A^— 
1, c(N) = A^+l}, for example. There exists a term Z\+Zn + \ since c(l) = l,c(N) = N+l. 
However, this term is equal to zero since it is one of the cases of Z^ + Z^+k = (70). 
Not only the excited states of Type I but also for a large number of low-lying excited 
states, Type III, for example, have terms Z^ + Z^ + k which eventually become zero. These 
states do not make any contributions to the relaxation dynamics. Instead, we find that 
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Table 2: Table of the asymptotic amplitudes vs total number of sites (alternating initial 
condition). 



i,k) 


\n[—B(EFP(i, k))] 


1,2) 


0.23825 - 


0.998191nM 


1,3) 


0.63787 - 


0.997361nM 


1,4) 


0.60321 - 


0.992481nM 


1,5) 


0.35850 - 


0.984021nM 


1,6) 


-0.01838 


- 0.971791nAf 


1,7) 


-0.49028 


- 0.956021nAf 


1,8) 


-1.03631 


- 0.936781nAf 


1,9) 


-1.64345 


- 0.914171nM 


1,10) 


-2.30316 


- 0.888251nM 



the second lowest excited state (Type II) determines the relaxation time r a i t for the case 
of the alternating initial condition. Note that the eigenvalue corresponding to the second 
lowest excited state obeys the KPZ scaling r a i t = — Re(l/7\4 2 nd) cx M 3 / 2 (see Section 4.2). 
Thus the local current asymptotically behaves as 

(h)t -> ^^fZY) + ^ndO^e^ 2 ^' as t -> oo. (84) 

The fitting from M = 256, 512, 1024 shows ln[£ 2nd (j fc )] = 0.23825 - 0.998191nM for both 
k odd (initially occupied) and even (initially empty). From this, one concludes £?2nd(iifc) oc 
M , which is the same with the step initial condition. 

Next, we analyze the EFP. Again, we find the lowest and third lowest excited states 
do not contribute. 

EFP{h k)t MKM-^^! 1 + B ( EFP h ^)> M2ndt as oo. (85) 

One observes B(EFP(i, k)) cx M" & , and @ k becomes smaller as the length k becomes 
longer. This behavior is similar with the step initial condition. 



7 Conclusion 

In this paper, we studied the long time asymptotics of the relaxation dynamics of the to- 
tally asymmetric simple exclusion process. We examined the local currents by the algebraic 
Bethe ansatz method. By evaluating the asymptotic amplitudes, we find the relaxation 
times starting from the step and alternating initial conditions are governed by different 
eigenvalues of the Markov matrix. The relaxation time of the step initial condition T s t ep is 
given by the nonzero eigenvalue of the Markov matrix with the first largest real part A4i s t 
as T s tep = — Re(l/A / list)- On the other hand, the relaxation time of the alternating initial 
condition r a it is given by the nonzero eigenvalue of the Markov matrix with the second 
largest real part A^nd as T a \ t = — Re(l/A42nd) for large number of total sites. The dif- 
ference between the step and alternating initial conditions is observed in another context: 
the current fluctuation. The GUE Tracy- Widom distribution appears for the case of the 
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step initial condition [17]. On the other hand, the current fluctuation for the alternating 
initial condition is described by the GOE Tracy- Widom distribution [18, 19]. Our results 
of the relaxation times is another aspect of the difference between the step and alternating 
initial conditions. 

In this paper, we examined the dynamics of the TASEP by use of the algebraic Bethe 
ansatz method. It is interesting to study other correlation functions and extend to other 
cases like open boundary, for example. The recent advances [44, 45] might be helpful for 
this case. It would also be valuable to study the link with the random matrix theory. The 
case examined in this paper is when the time is large enough, while the case obtained by 
the random matrix theory is when the time and the size of the system are comparable. 
These two results are considered to be supplementary to each other, and unifying the 
results by the Bethe ansatz would be an interesting problem. 
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